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A model for direct two band excitons in a disordered quantum well coupled to light in a cavity 
is investigated. In the limit in which the exciton density is high, we assess the impact of weak 
'pair-breaking' disorder on the feasibility of condensation of cavity polaritons. The mean-field phase 
diagram shows a 'lower density' region, where the condensate is dominated by electronic excitations 
and where disorder tends to close the condensate and quench coherence. Increasing the density 
of excitations in the system, partially due to the screening of Coulomb interaction, the excitations 
contributing to the condensate become mainly photon-like and coherence is reestablished for any 
value of disorder. In contrast, in the photon dominated region of the phase diagram, the energy gap 
of the quasi-particle spectrum still closes when the disorder strength is increased. Above mean-field, 
thermal, quantum and fiuctuations induced by disorder are considered and the spectrum of the 
collective excitations is evaluated. In particular, it is shown that the angle resolved photon intensity 
exhibits an abrupt change in its behaviour, going from the condensed to the non-condensed region. 

PACS numbers; 78.67.-n, 71.35.Lk, 71.36.+C, 42.50.Fx 



I. INTRODUCTION 

When the interaction between light and matter is 
strong, photons propagating in a dielectric can couple 
to excitons and form composite Bose particles known as 
polaritons 0. Since these composite bosons have a light 
mass, there has long been interest in the possibility that 
they might undergo Bose condensation at relatively high 
transition temperatures. While the concept of a Bose 
condensate of bulk polaritons has been discussed exten- 
sively 0i 13 ) in the open system, low-energy polaritons 
are merely long wavelength photons, which are not con- 
served 0- As a result, polaritons are unable to condense 
into the ground state making the bulk polariton conden- 
sate an intrinsically non-equilibrium phenomenon. How- 
ever, providing the lifetime of the polaritons is long com- 
pared to the thermalization time, the confinement of the 
photons in a microcavity may allow for the development 
of a quasi-equilibrium condensate. 

In recent years, improvements in the technology of 
semiconductor quantum wells has made the study of 
high-Q strongly-coupled planar microcavities almost rou- 
tine for III-V, II- VI and some organic semiconductors |^ 
Pumped both near resonance [1| (i.e. the cavity is 
excited at a 'critical' angle) or non-resonantly [fl, 13, llfl 
(i.e. the system is excited at large angles), sharp super- 
linear increase of the low energy polariton photolumines- 
cence has been observed. 

Polaritons are characterised by a very short radiative 
lifetime (of the order of ricak ~ 3 ps) due to leakage from 
the cavity. This has to be compared with the time scale 
to establish equilibrium with the lattice, riatticc, and the 
equilibration time of the polaritons themselves, Tpoiariton; 
the former controlled largely by exciton-phonon scatter- 
ing and the latter by exciton-exciton scattering processes. 
Dynamical (or quasi-equilibrium) condensation of polari- 
tons can be reached when Tpoiariton < Ticak < Tiattico- 



Due to the 'bottleneck effect' 0|, thermalization of small 
momenta polaritons by acoustic phonons is suppressed, 
causing a slow relaxation of polaritons in the low energy 
state. This process is found to be strongly suppressed at 
higher values of the pump power where superlinear emis- 
sion occurs 12J. Recently, evidence of quasi-equilibrium 
condensation, where the polariton life time reaches the 
relaxation time, including second order coherence in the 
optical field, has been observed . 

Within a bosonic picture involving tightly bound ex- 
citons, much work has been published on polariton dy- 
namics under non-resonant excitation, e.g. by making 
use of a rate equation approach In this paper, we 
will focus on the signatures of the fermionic nature of po- 
laritons when the excitation density increases. Here, the 
theoretical framework is not fully established even in the 
case of thermal equilibrium. The relative role of Coulomb 
forces (i.e. direct interaction between electron-hole pairs) 
and photon mediated interactions has not been studied. 
There is also an incomplete understanding of the effects 
of disorder, inelastic scattering and 'decoherence' phe- 
nomena, this being very important in drawing a distinc- 
tion between polariton condensation and conventional 
lasing. Furthermore, previous theories, including those 
that incorporate effects due to disorder and decoherence, 
have been at the mean field level. Our aim here is to 
address all three of these issues. 



A. Background 

Beginning with the pioneering work of Keldysh and 
coworkers 14j, the continuous transition between the 
Bose condensed and BCS-like phase of the equilibrium 
electron-hole system has been discussed within the frame- 
work of a mean-field theory characterised by an order pa- 
rameter involving the coherent polarisation. Refinements 
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of the mean-field theory by Comte and Nozieres |3 to 
include the effects of screening provided a consistent the- 
ory of the electron-hole liquid and condensate phases. In 
recent years, considerable efforts have been made to ex- 
plore the influence of light-matter interaction in the non- 
equilibrium electron system driven by an external phase 
coherent laser source — the optical Stark effect [l6|. 
However, only recently has the effect of quantum phase 
coherent coupling of photons and excitons in the closed 
system been addressed. 

Cavity polaritons interact both directly via the 
Coulomb coupling of their excitonic part, and indirectly 
via occupancy constraints on the excitonic component, 
sometimes referred to as phase-space filling. Since polari- 
ton masses are typically very small, the coherence tem- 
perature increases rapidly with number density. Once 
the characteristic interaction energy, either the dipole 
coupling energy or the Coulomb interaction, is larger 
than the coherence temperature of free polaritons, one 
must use an interacting theory. If the density is still 
low enough that excitons are not strongly overlapping, 
an appropriate approximation is to treat the excitons as 
bosons with a short range repulsion, or equivalently as 
two level systems. However, if the exciton density is high 
enough that the excitons become unbound then the sepa- 
rate electron and hole (fermionic) degrees of freedom are 
recovered. 

Note then that there are at least two and generally 
three crossovers as a function of increasing excitation: 
from a dilute gas of weakly interacting polaritons to a 
strong coupling regime, where nevertheless excitons are 
bound; then to a regime where the density of excitons is 
high enough that they unbind (and which may be weakly 
coupled to light, in the sense that the Rabi frequency is 
small in comparison to the kinetic energy). If, as is usual, 
the latter regime is one of weak coupling, nevertheless at 
high enough excitation one may reach a strong coupling 
regime where the optical field provides the largest scale in 
the problem. It is in the intermediate, fermionic, regime 
of weak coupling, that conventional semiconductor lasers 
operate [l^. 

Previous theoretical work has addressed the question 
of the existence of coherent, condensed polariton states 
in a simplified model. Taking just a single cavity photon 
mode, polariton condensation has been explored within 
the framework of an effective 'zero-dimensional' model in 
which excitons are trapped in isolated quantum dots |18| . 
Treated as two-level systems, the Hamiltonian for the 
coupled exciton-photon system can be presented as a 
Dicke model. The model allows for an inhomogeneous 
broadening of the levels, while a quasi- equilibrium con- 
dition was imposed by fixing the total number of excita- 
tions. At the level of mean-field (which becomes exact 
when the number of exciton states coupled by the same 
microcavity field is large), the ground state of the polari- 
ton system can be written as a variational wave function. 
This wave function is a superposition of coherent states 
of the exciton and photon in which both the amplitude 



of the cavity field and the electron polarisation acquire a 
non-zero expectation value. 

Within the same two-level system model, the effect of 
decoherence processes have recently been studied in . 
Treated again at the mean-field level, decoherence grad- 
ually suppresses the excitonic component of the order 
parameter, driving the system toward a gapless, weak- 
coupling regime similar to that of a conventional laser. 

Although the zero-dimensional system does provide in- 
tuition about the collective properties of the polariton 
condensate, by its construction, it is not able to discrim- 
inate between the BEC and BCS-like character of the 
condensate phase in the low and high density limits. Fur- 
thermore, the restriction of the model to a single cavity 
photon mode does not allow for the consideration of fiuc- 
tuations and dispersion of the excitations in the polariton 
condensate phase. Finally, in the two-level formulation, 
Coulomb interaction is not considered self-consistently, 
and both at mean-field level and beyond mean-field, the 
coherent polarisation and the coherent photon field are 
locked together. Among other consequences, this means 
that in this model the screening of Coulomb interaction 
by an increase of the density of excitations cannot be 
taken into account. On the other hand, such an effect 
is mimicked by the fermionic space filling effect, which 
also has the result of increasing the photon component 
at high excitation levels. 

In the following, we will tunc the strength of the 
light-matter interaction so as to explore the interplay 
of exciton and polariton condensation in a quasi two- 
dimensional microcavity system. Here, in contrast to the 
localised exciton theory, one can expect weak potential 
disorder to impact strongly on the integrity of the con- 
densate phase j_20j. Therefore, to keep our discussion 
general, we will include the symmetry-breaking effect of 
a weak potential disorder on the condensate. Although 
the mean-field theory has the capacity to describe the 
entire continuous interpolation between the low density 
BEC phase and the high density BCS-like phase, for sim- 
plicity, we will limit our analysis of the mean-field theory 
and fiuctuations to the high density phase. 



B. Model 

The many-body Hamiltonian for the coupled electron- 
hole/photon system can be separated into constituent 
components according to 

H — /l/Vox = Wei + Wdis + Wph + 'Hint ■ (1) 

Here, TLd represents the interacting Hamiltonian of the 
electron and hole degrees of freedom in a semi-conducting 
quantum well, while Tidis describes the random poten- 
tial generated by quenched impurities and lattice defects. 
Tiph denotes the Hamiltonian associated with the cavity 
photons and finally TL-mt represents the exciton-photon 
dipole coupling interaction. 
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Following the notation introduced in Ref. [23 
adopting an 'electron picture' in which bp and &t 



and 



and Oj^) represent the annihilation and creation opera- 
tors of an electron with momentum p in the conductance 
(valence) band (see Fig. the Hamiltonian for a direct 
band-gap semiconductor is given by 

^ci-5]ep(^^^^ + ap4) 
p 

q^^O \ p p / 



where 



qbp - 



P"p+q 



) is the total electron 



Pq = Epv _ _ 

density operator. Here, for simplicity, we suppose that 
the electrons are spinless. Taking the middle of the gap as 
the energy reference, we assume a parabolic conduction 
and valence band of particles 
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where, for simplicity, we have assumed the electrons and 
holes exhibit the same effective mass, m. At the level of 
mean-field, one can confirm that this assumption is in- 
nocuous and corresponds to replacing the effective mass 
with the double of the reduced mass. In contrast, devia- 
tion from parabolicity in the physical system impairs the 
capacity for perfect nesting of the Fermi surface and re- 
duces the potential for pair correlation. Finally, for con- 
venience, we have introduced in Eq. Q the concept of an 
effective Fermi energy sp = p^p/2m = {^j, — Eg)/2, which, 
along with the chemical potential fi, is determined self- 
consistently as a non-trivial function of the total number 
of excitations in the system (jSJ (see later in section llil A|l . 

In addition to the lattice potential, the electrons and 
holes experience a Coulomb interaction. In the high den- 
sity regime peio^ 3> 1, where pei — pj^/An denotes the 
areal density of the particles, (i.e. when the the excitonic 
Bohr radius oq = 2eo/e'^m is much larger than the av- 
erage distance between electron and holes the 
Coulomb interaction is screened due to both electrons 
and holes [Uli^: 



eo(q2 



In the two-dimensional system, the screening length is 
approximately set by the Bohr radius 1/k = cq. Hence, 
in this limit, the Coulomb interaction can be replaced by 
a short range contact interaction j23|, 



p' 



^ I y \pf -p\ < K 

\p'\=pf I \pf - p\ > K 



with a strength given by the angular average over the 
Fermi surface 



^ = 7^ya^ 
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FIG. 1: Schematic picture of the valence (a) and conductance 
(6) band, and the interactions included in the model Q. 



Here TZy — e^/2eoao denotes the Rydberg energy associ- 
ated with the exciton and sets the characteristic energy 
scale of the interaction. Therefore, as a consequence of 
the screening, an increase in the total number of elec- 
tronic excitations leads to a reduction in the effective 
strength of the Coulomb interaction. Notice that, al- 
though the condensed state may have a gap, we shall be 
discussing the dense limit, where the coherence length is 
much longer than and so this will not influence the 
important short range part of the potential. 

To incorporate the effect of disorder, we suppose that 
the electron-hole system is subject to a generic weak im- 
purity potential, 

Tidis- j dv [Ub{vp{v)b{r) + Ua{v)a\v)a{v)\ . 

Naively, since the electrons and holes carry opposite 
charge, their response to potential impurities should be 
equal in magnitude and opposite in sign. However, dif- 
ferences in the dielectric properties of the conduction 
and valence band electrons, combined with the potential 
asymmetry imposed by the quasi two-dimensional geom- 
etry, implies that the potentials Ua and Ub are not per- 
fectly anticorrelated. Furthermore, there are mechanisms 
of disorder such as variations in the width of the quantum 
well through atomic terracing, strain fields, etc., that act 
on electrons and holes with the same sign. Therefore, in 
the following, we will suppose that a generic disorder po- 
tential is comprised of two channels: A potential ?7„(r) = 
[J7f,(r) — Ua{v)]/2 which acts symmetrically as a 'charge- 
neutral' potential and, as we will see, does not effect the 
integrity of the condensate phase. A second, statistically 
independent, potential J7c(r) = [^/^(r) -I- C7a(r)]/2 acts 
as an asymmetric 'charged' component, and presents a 
pair-breaking perturbation which gradually destroys the 
excitonic condensate. As mentioned above, its effect in 
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the dense exciton system on the integrity of the exciton 
insulator phase has been explored in an early work by 
Zittartz [23. In both cases, we will suppose that the 
potentials Un.c are drawn at random from a Gaussian 
white-noise distribution, with zero mean, and variance 
given by 

{UnA^l)UnA^2)) = ^(ri - Ta) , (4) 

where r„ and Tc are the associated scattering times, and 
u = m/27r is the two-dimensional density of states of 
the spinless electron system. Restricting attention to 
the high density phase, setting spTn^c ^ 1, we will fur- 
ther suppose that the impurity potential impose only a 
weak perturbation on the electron and hole degrees of 
freedom, allowing therefore the dynamics to be treated 
quasi-classically. This completes the description of the 
electron-hole system. 

Free photons in the cavity are described by the micro- 
scopic quasi two-dimensional Hamiltonian 

P 

where their dispersion, w(p) = ^Ju^ + (cp)^, is quan- 
tised in the direction perpendicular to the plane of the 
cavity mirrors, lOc — cnj/La^, where L^i denotes the dis- 
tance between the mirrors, and j € Z. In the follow- 
ing, we will suppose that the electron-hole system en- 
gages just a single sub-band for which ujc 7^ 0. In prac- 
tice, providing the sub-band separation is substantially 
larger than the quasi-particle energy gap that develops 
in the electron-hole system, neighboring sub-bands can 
be safely neglected. 

Finally, in the dipole or 'rotating-wave' approximation, 
the photons are assumed to be coupled to the electron- 
hole system through a local interaction, 

Hint =9jdr [VXr)fet(r)a(r) + h.c] . 

Here, terms which do not conserve the number of exci- 
tations (i.e. the ones describing spontaneous creation or 
annihilation of a photon and an exciton) have been ne- 
glected. On one hand, we expect this approximation to 
be valid in the weak coupling limit between photon and 
matter, where these additional terms can be 'integrated 
out', giving a small 'renormalisation' of the single parti- 
cle electronic and photonic dispersion laws. At the same 
time, since we are dealing with a system where the tem- 
perature is much smaller than the photon frequency, we 
can neglect the tiny spontaneous population that would 
be generated by these non-resonant terms. 

To mimic the effect of the external photon source, 
we suppose that the electron- hole/photon system is held 
in quasi-equilibrium by tuning the chemical potential /i 
in ^ to fix the total number of excitations 

iVex = 5I^P^P + ^II(^P^+«P°p) ■ 
P P 




FIG. 2: Relative number of (total) electrons rid = Pci/Pox and 
(condensed) photons Uph ~ Pph/pcx are plotted against pcxflo 
for fixed values of the coupling constants g — 1 and Ac ~ 10 
and for d — 200. Note that, for pcxflo > 25, the system enters 
a strong coupling regime, where our model cannot be applied. 
In the inset, the number of photons PphO-o is compared with 
the number of electrons in the condensate, A'^|e| ~ \'S\i'aQ (see 
Eq. I21ll 'l. For Poxflo > d/An ~ 16, the excitations in the 
condensate are dominated by photons. 



However, how the system chooses to portion the exci- 
tations between the electron-hole and photon degrees of 
freedom depends sensitively on the properties of the con- 
densate. 

This completes the formal construction of the micro- 
scopic many-body Hamiltonian. In order to explore the 
mean-field content and the collective excitations of the 
light-matter system, we will draw on the theory of the 
weakly disordered superconductor. Using conventional 
field theoretic methods, we will show that the low-energy 
properties of the system can be cast in the framework 
of an effective field theory with an action of non-linear 
sigma-model type (21| . Partitioning our discussion into 
properties of the mean-field content of the theory and the 
role of fluctuations, we will focus on the phase diagram of 
the system as a function of density and disorder on one 
side, and the nature of the collective excitations on the 
other. Since the field theory technology disguises much of 
the physical content, we will close the introduction with 
an overview of the findings of the mean-field theory. As 
well as summarising the main results, it will provide a 
useful prospective for the formal analysis. 

In the high density regime, the mean-field content of 
the action will be shown in section IIIII to be insensi- 
tive to the 'neutral' component of the disorder poten- 
tial, J7n(r), a manifestation of the Anderson theorem 
in the exciton context. This reproduces the mean-field 
thermodynamic properties of a system subject only to a 
'charged' weakly disorder potential, Uc{r). By analysing 
these mean-field equations, the generalisation of those 



derived by Zittartz |20j for the exciton insulator system, 
we will describe the rearrangement of the ground state 
due to condensate formation. In particular, we will show 
that the condensate of the polariton system is charac- 
terised by an order parameter which engages the combi- 
nation of the phase locked polarisation and photon ampli- 
tudes. The mean- field phase diagram (see Fig. O can be 
characterised in terms of the total density of excitations, 
PcxOo, the disorder strength ^ = l/rc7^y and the follow- 
ing dimensionless material parameters, respectively the 
coupling to the photon field and the Coulomb coupling 
strength (see expression (jSJ): 



(6) 
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,(x) = gcvL^ = 



y/l + 4{d-x) 



de- 



Here, the constant of proportionality Ac = 4ao/Lj 
pends on the thickness L± of the quantum well, x = 
{lUc — a*)/ Ti-y denotes the characteristic energy sepa- 
ration between the photon band edge and the chemi- 
cal potential (measured in units of the Rydberg), and 
d = {ujc — Eg)/ TZy > 0. Note that, as in BCS theory, gc 
scales with the volume, i.e. in the thermodynamic limit 
gcL^ is finite while, by contrast, g scales with the square 
root of the volume, i.e. g'^L'^ is finite. 

When the density PexOo l^^i ^^'^ majority of excita- 
tions are invested in electrons and holes (see figure|2J) and 
therefore, as the density increases, the chemical potential 
rises linearly. Here, providing the density is large enough 
to keep the particles unbound, the condensed phase is 
reminiscent of the exciton insulator one, even if a small 
fraction of photons do contribute to the condensate (see 
the inset of figure [2J|. Eventually, as the chemical po- 
tential approaches the band edge, lOc, the photons are 
brought into resonance and the character of the con- 
densate changes abruptly. Here, the excitations become 
increasingly photon-like, with PcxO-Q diverging exponen- 
tially when the chemical potential converges on lOc- At 
the same time, screening suppresses the excitonic cou- 
pling constant, determining the majority of the excita- 
tion in the condensate to acquire a photon-like character. 
This is clearly shown in figure [3 where the normalised 
number of total electrons Poi/Pox and photons Pph/pcx 
(see Eq. (I2UII ') are plotted as a function of the density 
of excitations. In the inset of the picture, instead, the 
fraction of photons to excitons in the condensate (for the 
exact definition, see later on, Eq. l(2U) clearly shows an 
abrupt growth as the density of excitations in increased 
above the value d/An, which we will see to coincide with 
the maximum density which can be reached at /i = cjc in 
absence of photons. 

Staying at the level of mean-field, the 'charged' compo- 
nent of the disorder potential presents a symmetry break- 
ing perturbation which depletes the excitonic component 
of the condensate. As such, its effect on the condensate 
depends sensitively on the density of excitations. When 
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FIG. 3: Mean-field phase diagram for g = 1, Xc — 10 and 
d = 200; pe^a^, represents the excitation density, while ^ = 
l/TcTly the disorder strength. 



1 < PcxCto d/47r, the condensate has a particle-hole 
character and the phase diagram (Fig. mimics the be- 
haviour of the symmetry broken exciton insulator |2fl| | - 
and therefore of a superconducting system in presence of 
magnetic impurities. When the scattering rate Tc is com- 
parable to the value of the unperturbed order parameter, 
the system enters a gapless phase before the condensate 
is extinguished altogether. At larger densities, the con- 
densate becomes photon-dominated and therefore robust 
against any value of the disorder potential. Interestingly, 
however, the residual effect of the disorder exposes a sub- 
stantial region of the phase diagram where the system is 
gapless. Here, the condensate acquires all of the conven- 
tional characteristics of a semiconductor laser — i.e. a 
substantial coherent optical field, but a gapless spectrum 
of electron-hole pairs with negligible electronic polarisa- 
tion l^^. Once again, at sufRciently large densities, the 
system enters a strong coupling regime, where the valid- 
ity of the model becomes doubtful. 

This completes our summary of the mean-field content 
of the theory. The remainder of the paper is organised as 
follows. In the next section, we develop a field theory of 
the weakly disordered polariton system. Focusing on the 
saddle-point approximation, in section IIIII we elaborate 
and expand our discussion of the mean-field content of 
the model. Using these results as a platform and expos- 
ing analogies with the symmetry broken superconducting 
system, in section Hvl the attention is then turned to the 
impact of thermal, quantum and disordered driven fluc- 
tuations. Here, we explore the spectrum of the collective 
excitations and, in particular, we discuss how the tran- 
sition to the condensed phase provide a signature in the 
photoluminescence. Finally, in section we conclude. 
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II. PATH INTEGRAL FORMULATION 

The development of a low-energy field theory of the 
weakly disordered exciton-light system mirrors closely 
the quasi-classical theory of the disordered superconduc- 
tor. By presenting the quantum partition function of the 
many-body system as a coherent state path integral, the 
ensemble or impurity averaged free energy can be cast in 
terms of a replicated theory. At the level of the mean-field 
theory, this approach is equivalent to the self-consistent 
Hartree-Fock treatment where the impact of disorder is 
treated in the self-consistent Born approximation. Using 
the mean-field result as a platform, the field theoretic 
formulation provides the means to develop a field theory 
of the low-energy excitations of the polariton condensate. 

The quantum partition function Z — tr[e^'^*^^~^^™)], 
where /3 — is the inverse temperature, can be ex- 

pressed as a coherent state path integral over fermionic 
and bosonic fields. To facilitative impurity averaging of 
the free energy over the disorder potentials Q , it is con- 
venient to engage the 'replica trick' p^ . 



F = 



4 lim ■ 

(3 n 



1 



Although the analytic continuation presents difficulties in 
treating non-perturbative structures, the technique can 
be applied safely to the mean-field and collective proper- 
ties discussed in this work. Once replicated, a Hubbard- 
Stratonovich decoupling of the Coulomb interaction ob- 
tains the field integral for the quantum partition function 



exp 



dr J dr (j)\r, r) (^dr + i?) 0(r, t) 



(7) 



where, in order to lighten the notation, wc have omitted 
the replica index carried by all of the fields. Here, 5ph 
and Ss correspond respectively to the bare photonic and 
excitonic components of the action, 

5ph= / drJ2rpir)[dr+Lo{p)-fAMr) (8) 
^" P 

^s-^^ drj dr |E(r,r)p , (9) 



while, arranging the fermionic fields into a Nambu-likc 
spinor (j)^ = {b^,a'^) and cj) = (&, a)^, the (single-particle) 
action for the internal electronic degrees of freedom 
(which, with an abuse of terminology, we will later on re- 
fer to as 'particle-hole' components) takes the form rem- 
iniscent of the Gor'kov or Bogoliubov de-Gennes super- 
conducting Hamiltonian, 



Here, ai represent Pauli matrices operating in the 
particle-hole subspace, and A = |A|e*'' represents the 
complex 'composite order parameter' 







lAkie- 



H 



ip + Un(r) a3 + Ucir) + A 



(10) 



involving both the polarisation E and photon tp fields. 
Taking into account the contact nature of the screened 
Coulomb interaction, the coupling constant gc can be 
derived starting from the expression iPJ and is given in 
Eq. (|HJl. At the level of mean-field, the excitonic order pa- 
rameter signals the development of the anomalous pairing 
amplitude 

E{r) = -g,L^g.s.\b\r)air)\g.s.) , 

while the ordinary Hartree-Fock pairings, 
(g.s.|6'''(r)6(r)|g.s.) and (g.s.|a(r)a'''(r)|g.s.), crucial 
in the low-density limit, effect only a small renormal- 
isation of the single-particle energy ^p. Restricting 
attention to the high density regime, such renormalisa- 
tions can be neglected (or absorbed into a redefinition 
of the position of the band edge). 

As noted above, the single-particle Hamiltonian for 
the electron and hole mirrors the structure of the quasi- 
particle Hamiltonian for a superconductor. In particular, 
in the absence of disorder, an electron of energy in the 
valance band can pair with a hole of energy ^p in the con- 
duction band. The development of a homogeneous order 
parameter A induces a gap in the quasi-particle density 
of states at the effective Fermi level ep- Already at this 
level, the insensitivity of the condensate to 'neutral' dis- 
order Un is apparent: Even in the presence of the disor- 
der potential, the electron with a single-particle energy 
can pair with a partner hole with the same energy and 
form a homogeneous condensate, a reflection of the An- 
derson theorem of superconductivity in the electron-hole 
system. By contrast, the symmetry breaking potential 
Uc lifts the degeneracy of the electron and hole degrees 
of freedom and imposes a pair-breaking perturbation on 
the condensate. As a result, one expects a gradual de- 
pletion of the condensate. In fact, we will see later on in 
section UlI Bl when establishing the phase diagram, that, 
since photons are almost insensitive to the effect of dis- 
order, only the 'excitonic contribution' to the condensate 
will be affected, while at very high densities, the photonic 
component will be able to restore coherence. 

In principle, one could explore the mean-field content 
of the coupled system by developing a diagrammatic self- 
consistent Hartree-Fock approximation. However, in the 
following, we will be concerned with the impact of fiuc- 
tuations on the integrity of the condensate phase. As 
well as thermal and quantum fluctuations, mechanisms 
of quantum interference affect the spectrum of quasi- 
particle excitations above the condensate. To properly 
accommodate these effects, it is helpful to construct an 
effective low-energy theory which encodes the diffusive 
character of the dynamics. To classify the different chan- 
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nels of interference, it is helpful to identify the symmetry 
properties of the coupled system. 

As well as its invariance under the global gauge trans- 
formation 



P*x''^3/2 



-ix' 0-3/2 



■0 1-^ ipe^^ 



(11) 



in the absence of the symmetry breaking potential C/c, the 
matrix Hamiltonian H1U|I exhibits the discrete particle- 
hole symmetry H — — CT2-ff^0'2- Its presence facilitates 
mechanisms of quantum interference which effect the low- 
energy properties of the quasi-particles excitations. Fol- 
lowing a standard procedure (see, e.g., Ref. 24] ), in or- 
der to properly take into account for the 'soft diffusion 
modes' associated with this symmetry, it is useful to dou- 
ble the field space setting 



* V2 \'72 '"^ 



V2 



0-2 



As a result, the fields $ engage a total of 4 x n compo- 
nents (i.e. n replica, particle-hole and 'charge conjuga- 
tion' cc), while the redundancy implied by the field dou- 
bling imposes the symmetry relations $ = (727$^ and 
$ — —$^(727^, with 7 = It is important to note 

that the impact of the full charge-conjugation structure 
of the fields is manifest only in the mesoscopic fluctua- 
tion phenomena which effect only the low-energy quasi- 
particle excitations. In particular, such mechanisms of 
quantum interference do not effect the properties of the 
system at the level of the saddle-point or mean-field. 

Cast in the form {Tj), the action mirrors closely that 
studied in the context of the disordered superconduc- 
tor and, to identify the effective low-energy theory, one 
may draw on the existing literature. In particular, in the 
quasi-classical and dirty limits. 



1 1 ,A, 

e_F > — > — , |A| 



(12) 



it has been shown [2J, [23 that the long-range properties 
of the system can be expressed in terms of an action of 
nonlinear tr-model type. 



(Z") = J DiiP^ij) J £>(!]*, S)e-^p"~^^ J DQ, 



where 



S[Q] = ^ j drtv[D{^Qf 



-4 ( €^3^^ - i Aaa ]Q-- (fTgaf Q)' 



(13) 



T are compatible with the 'charge-conjugation 'symme- 
try properties which Q inherits from the dyadic product 
(T3<i> (g) I*. Qualitatively, the fluctuations of the matrix 
field Q encode the effect of long-range mechanisms of 
quantum interference. 

Combined with the bare actions for the auxiliary 
fields © and ©, the non-linear cr-model action fully 
encodes the effect of weak disorder in the cavity polari- 
ton system. In the absence of the order parameter, the 
action describes the diffusive density relaxation of, sep- 
arately, the electron and hole degrees of freedom. The 
development of the order parameter |A| describes the re- 
arrangement of the ground state due to the formation 
of the polariton condensate. Finally, scattering in the 
channel of the 'charged' disorder potential projects out 
the diffusion modes which act coherently on the particle- 
hole degrees of freedom. In order to explore the capacity 
of the system to undergo condensation, it is necessary 
to investigate how the order parameter, and symmetry 
breaking impurity potential revise the structure of the 
saddle-point solution. 



III. SADDLE-POINT ANALYSIS: MEAN-FIELD 
THEORY 

The mean-field content of the theory is encoded in the 
saddle-point structure of the effective action. Varying 
the total action (|13|l with respect to subject to the 
nonlinear constraint = 1, the saddle-point equation 
takes the form 

£»V(QVQ) + [Q, eVgfrf - i Acts] 



2r, 



■[Q,cj^afQa:,af]^Q. (14) 



Such equations of motion are familiar from the quasi- 
classical theory of the weakly disordered superconduc- 
tor, where the field Q is identified as the average quasi- 
classical Green function |26[. The second term describes 
the 'rotation' of the average quasi-classical Green func- 
tion due to the formation of a condensate A, while the 
last term describes the pair-breaking effect of the pertur- 
bation imposed by the symmetry breaking disorder. 

Similarly, varying the action with respect to the exci- 
tonic order parameter field S, and the photonic field "^j 
one obtains the coupled self-consistency equations: 



9, 



(r) 2^ tr [(T_Q„+h .„(r)] 



4/3 



i^e'P""[iw,,+w(p)-/i]V'c..,p , (15) 



Here, D — VpTn/2 represents the diffusion constant asso- 
ciated with the 'neutral' disorder, while the matrix field 
Q(r) = r(r)A(T3(T3'^T~^(r) carries replica, time and in- 
ternal (particle-hole and CC) indices, and the generators 



where the matrix (T_ — ai — iai projects onto the off- 
diagonal particle- hole channel and where ujh — 2Z7r//3 
are the bosonic Matsubara frequencies. Since the pho- 
tonic and excitonic interactions are decoupled in the same 
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FIG. 4: The order parameter |A|/|Ao| (solid line) and energy 
gap egap/|Ao| (dashed line) versus the strength of disorder 
1/tc|Ao|. Here, |Ao| represents the value of the order pa- 
rameter in absence of disorder. The critical value at which 
the order parameter disappear is given by l/rc|Ao| = 1/2, 
while the gapless region starts at l/rc|Ao| — e"'^''* ~ 0.46 (or 
equivalently at ( = 1/tc\A\ — 1). 



channel, it is not surprising that, at the level of the 
saddle-point, the two order parameters cannot be inde- 
pendent but instead are coupled together by the con- 
straint implied by (|15|l . This means that in the con- 
densed phase both order parameters develop in concert. 
In other words, in the condensed phase of the present mi- 
croscopic theory, both electrons and photons participate 
in the formation of the condensate — it is only their rel- 
ative proportion that depends (sensitively) on the total 
number of excitations iVox present in the system. 

The mean-field equations H14|) and (|15|l can be solved 
simultaneously with the mean-field Ansatz that the two 
order parameters are space and time independent (i.e., 
E(r,T) — T, and ipuj^.p = V"5a)h,oi5p,o)j and the matrix 
field Q{r) is homogeneous in space and diagonal in the 
Matsubara indices. Now, through the transformation 



Q 



o2XCT3/2r)„-ix<^3/2 



, the (homogeneous) global phase 



X of the composite order parameter A can be gauged out 
of the action H13|l and the equation of motion (|14() can 
be solved, giving 



0-3(73 



(72 



where the variable u satisfies the equation [22 

En _ 1 u 

lAI 



(16) 



As expected, in the mean-field approximation, the so- 
lution of the quasi-classical equations of motion become 
insensitive to the symmetry preserving disorder potential 
(i.e. the solution is independent of the classical diffusion 
constant D). Such independence is a manifestation of the 



Anderson theorem that protects the integrity of the con- 
densate to impurities which do not lift the particle-hole 
symmetry of the system. 

In the homogeneous case, the self-consistency equa- 
tions for the order parameters 1)15(1 require the phases of 
the excitonic order parameter S and of the photonic field 
■0 to coincide, so that |A| = |E| -I- g\'i\}\ and from which 
follows the constraint: 



1. 
9 



9c 



(17) 



Indeed, the same restriction has been found in Refs. 0, 
19] in the context of the Dicke model for the localised 
exciton-photon system. While the chemical potential fi 
does not exceed the cavity mode edge lUc, an order pa- 
rameter |A| ^ is developed, to which both photons and 
the coherent polarisation contribute. Once cast in terms 
of the effective order parameter A, the remaining saddle- 
point or mean-field equations assume the canonical form 
first reported by Abrikosov and Gor'kov in the context 
of the symmetry broken disordered superconductor |27j| . 
and later, in the present context, applied to the symme- 
try broken exciton insulator |20j |. To make this corre- 
spondence explicit, one may define the composite order 
parameter through, say, the photon condensate fraction, 
and an effective coupling constant as 



1 + ^{ujc - m) 



9 



9\M 



9cs = gc 
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after which the self-consistency equation reads 
lAI 



TT X - 1 

« 2-^ .FTT 



(18) 



In the absence of symmetry-breaking disorder, u = 
e„/|A| and the self-consistency equation coincides with 
the 'gap equation' for a bulk BCS superconductor. Keep- 
ing this analogy, convergence of the Matsubara summa- 
tion demands the inclusion of an energy cut-off Ec as- 
sociated with the contact Coulomb interaction, analo- 
gous to the Debye frequency scale in the superconduc- 
tor, which sets the overall scale of the order parameter. 
In the present case, the momentum exchange associated 
with the Coulomb interaction is limited by the screening 
length K. As a result, the contact interaction, and there- 
fore, the Matsubara summation, must be cut-off at the 
energy scale 



Er 



ppK 



PF 



(19) 



In the high excitation limit where the photon coupling 
dominates, the cutoff should be substituted with ep, the 
overall bandwidth for electron-hole pairs. 

At zero temperature, an explicit solution of the gap 
equation reveals the gradual suppression of the conden- 
sate as a function of the disorder strength 1/tc|Ao|, lead- 
ing to a complete suppression of the order parameter 



n 



when 2/tc = |Ao|, where |Ao| denotes the order pa- 
rameter m the unperturbed system (see Fig. 0)). (For 
a general review in the context of the symmetry bro- 
ken superconductors, see, e.g., Ref. (2^.) As with the 
disordered superconductor, the quasi-particle energy gap 
goes to zero faster than the order parameter, exposing 
a gapless phase when 2/tc > e~'^/^|Ao| or equivalently 
when the dimensionless parameter <^ = 1/tc|A| exceeds 
unity. To understand how the phase behaviour of the 
Abrikosov-Gor'kov equations translate to the polariton 
system, it is necessary to explore the dependence of the 
chemical potential on the number of excitations. 



A. Chemical Potential and Number of Excitations 




Hitherto, we have absorbed the chemical potential 
into the definition of an effective Fermi energy, sf — 
{fx — Eg)/2, assuming this to be the largest energy scale 
in the problem. However, the chemical potential /i is 
itself fixed by the total number of excitations in the sys- 
tem In the following, we will use this condition to 
explore the parameter range of validity of the approxima- 
tions used to construct the theory, and to investigate the 
phase diagram, treating TVox and the symmetry breaking 
disorder potential I/tc as independent variables. 

The validity of the model, and therefore of the self- 
consistency equation (|18|l . is restricted by the hierar- 
chy of energy scales contained in Eq. (|12|l . Moreover, 
as implied by Eq. H19|l above, the energy cut-off of the 
Coulomb interaction is, itself, dependent on the chemical 
potential according to the relation = y/TZy{iJ, — Eg), 
where TZy is the Rydberg energy of the electron-hole sys- 
tem. Taken together, the validity of the mean-field the- 
ory requires ep = — Eg)/2 > Ec (i.e. (/z — Eg)/ TZy > 
4), while to stay in the weak coupling limit, we must have 
Ec ^ |A|. In this limit, at zero temperature, the total 
number of excitations iVox, which is proportional to the 
partial derivative of the free energy respect to the chem- 
ical potential, can be expressed in terms of the chemical 
potential as it follows: 



FIG. 5: Normalised density of excitations pcxOo ^ function 
of the rescaled chemical potential x = {uc — fJ^)/ TZy (in loga- 
rithmic scale) for d = {luc - Eg) /TZy = 200, g = 1, Xc = W 
(bold line) and for g = (dashed line). The value of the 
density at a; = for g = 0, d/A-K ~ 16, has been marked. In 
this limit, all excitations appear as electrons and holes, al- 
lowing the population of states with energies in excess of loc- 
The dependence from of the dimensionless effective coupling 
constant uL^gaa on pcxtso in both cases is shown in the in- 
set. Note that the weak coupling limit imposes the constraint 
Pcxflo < 25. 



ep > Ec^ |Ao| obtains the BCS solution 



|Ao| 



Er 



sinh (l/vL'^gcs) 



Staying within the weak coupling regime, vL^g^s < 1, 
to check that the system remains in the high density 
phase, let us recall that the Coulomb coupling constant 
is a function of the chemical potential (see Eq. Q). It is 
convenient to introduce the dimensionless coupling con- 
stants (|nj. With these definitions, the total density of 
excitations PcxOo = ^c^l Li^a^, takes the form 



212 I 



vL 



2/" 



En 



The first term describes the total number of photons 
■^ph = while the second represents the total number 
of particle-hole excitations in the system, Mc\ ~ I'L^ep- 
Here we have neglected terms of order 0{\A\/eF) associ- 
ated with the rearrangement of states close to the Fermi 
energy due to the formation of the condensate. Taking 
into account the non-algebraic dependence of the order 
parameter as well as the coupling constant gc on /i, the 
variation of the chemical potential with A^ox can be in- 
ferred only numerically. 

For simplicity, let us first consider the case of a clean 
quantum well, I/tc = 0. In this solution of the 

zero temperature self-consistency equation in the limit 
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where = Afph/L'^ denotes the photon density and 
Pel — A/'ci/L^ = l/r^ denotes the total electron- hole den- 
sity (and from which one can infer the value of the 
electron- hole system). Inverted, Eq. H20|l describes the 
variation of the chemical potential (through x) as a func- 
tion of the material parameters g, Ac, d and the exciton 
density pox- 

In figure \S\ the (normalised) density of excitations 
PoxIq plotted against the renormalised chemical po- 
tential X. When the density is low, the vast majority of 



excitations are invested in particle-hole excitations and 
the chemical potential scales linearly with the density, 
reflecting the constant density of states in the quasi two- 
dimensional system. As the density grows, the chemical 
potential /i converges on the photon band edge o^c, i.e. 
the variable x diminishes to zero. Here, close to the res- 
onance, the excitations become increasingly photon-like 
and Po-^Uq diverges exponentially. The maximum density 
of electron-hole excitations which can be added at /i = Wc 
in the absence of the photon interaction, g — 0, is given 
by d/An. This value plays an important role, since, it 
represents the value of the density at which the conden- 
sate becomes dominated by photons. In fact, although 
both photons and electrons contribute to the formation 
of the condensate (i.e. to |Ao|), the electronic contri- 
bution is dominant at low densities, while the screening 
of the Coulomb interaction results in the growth of the 
photon contribution at large densities. Such behaviour is 
illustrated most clearly in the inset of Fig. [S] where the ef- 
fective coupling constant vLF'geS is plotted against Pcxfto- 
Increasing the density of excitations above d/Aw, the ef- 
fective coupling constant abruptly starts to increase until 
one eventually reaches the strong coupling regime where 
the validity of the microscopic Hamiltonian becomes un- 
certain. These characteristics are reflected in flgure 
where the fraction of total electrons Pci/Pcx and photons 
Pph/Pox are plotted as a function of the density of exci- 
tations. Moreover, making use of the equality (|17|l . the 
inset shows the variation of the fraction of photons to 
excitons in the condensate, 

yPpha'p ^ g_V^ ^21) 

where = |E|t/aQ. Note that this ratio starts to devi- 
ate from zero as soon as particles are added to the system. 
In other words, even if small, there is a contribution from 
the photons to the condensate at 'low densities', when 
PcxO'Q ^ d/Air. Moreover, as already noted in the intro- 
duction, the increase of this quantity with the increase 
of the total density of excitations is due, on one hand, 
to the increase of the photon content as one approaches 
the resonance while, on the other hand, is related to the 
gradual decrease of Coulomb interaction and increase of 
the exciton-photon interaction. 

So far we have focused on the character of the polari- 
ton condensate. In the next section, staying at the mean- 
field level, we will explore its integrity when effected by 
the symmetry-breaking disorder potential. Note that we 
have chosen particular values for the parameters d, Ac 
and g in such a way as to expose a weak coupling regime 
in the high density limit. Changing their values, the main 
characteristics of the phase diagram remain valid, but the 
specific numbers can be varied. In particular, changing d, 
the regions where electron-like (lower density) or photon- 
like (higher density) excitations are dominant can be en- 
larged or reduced in size. The range of the condensed 
region for pexOo < d/An can be varied through Ac, while 
changes in g characterise, when Poxio ^ d/Air, variations 
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FIG. 6; Order parameter jA|/7^y versus p^^a^S for g = 1, 
Ac = 10 and d — 200 and different values of the disorder 
strength ^ = I/tcTZj. 

in the gapless region and changes in the value at which 
the strong coupling regime is entered. 



B. Phase Diagram 

In the symmetry broken disordered system, the pho- 
ton density (as recorded in the magnitude of the order 
parameter) has to be rescaled by the weight |A|/|Ao|. 
Fig.0]shows the variation of the order parameter against 
the disorder strength l/rc|Ao|. Taking into account the 
dependence of the unperturbed order parameter |Ao| on 
the chemical potential (and therefore on the number of 
excitations), |A|/7?.y is plotted in figure El against pcxOo 
for different values of ^ = 1/TcTZy. 

Let us consider first the dependence of the order pa- 
rameter on the excitation density in the clean case: Start- 
ing from poxflo ^ 1' 1^0 1 first slightly decreases, due to 
the switching off of the Coulomb interaction, and then 
starts to increase rapidly as 

\Ao\/ny~V2^g^Jpexal , (22) 

for pcxO,Q > d/An. At these high values of the density, 
the excitations in the condensate are mainly photon-like. 
Therefore, increasing the disorder strength, it is not sur- 
prising that the profile of the order parameter changes lit- 
tle. In can be shown analytically that, in this limit, the 
condition 1/tc|Ao| > 1/2 cannot be reached. At lower 
densities, however, where the excitations are mainly 
electron-hole like, the order parameter starts to close un- 
til, increasing the value of disorder, it eventually van- 
ishes everywhere. This behaviour is more evident from 
the shape of the phase diagram (see Fig.|2Jl: In the region 
of density where the excitations are mainly electron-hole 
like, the condensed phase gets reduced, while, when the 
photon content starts to dominate, the condensed region 
expands until eventually, when PcxOo — d/Air, the disor- 
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FIG. 7: Three-dimensional plot of the critical temperature, 
Tc/lZj, versus the disorder strength, ^ = I/tcTZj, and the 
excitation density, pexflo id = 1, Ac = 10 and d — 200). The 
contour line where the critical temperature goes to zero sep- 
arates the zero-temperature condensed phase from the non- 
condensed one and it coincides with the contour given in fig- 
ure El 

der becomes ineffective, and the line separating the two 
phases becomes asymptotically straight. 

Analogously, one can evaluate the density and disorder 
dependence of the energy gap 

/ , \ 3/2 

e,.p = |A|(l-C^/^) 

and determine the region where the excitation spectrum 
becomes gapless, egap = (which can be equivalently es- 
tablished from the condition C, = 1/tc\A\ > 1). As a 
result, in the region dominated by the electronic excita- 
tions the gapless phase is very small and lies just below 
the non-condensed phase, while, when pcx^o ^ d/Air, the 
line which separates the condensed gapless region from 
the gapped one grows approximatively as V^g^/p- 

C. Finite Temperature 

Although a rigorous treatment of the finite tempera- 
ture mean-field theory would require an analysis of the 
thermal population of the photon degrees of freedom, 
one can gain a qualitative insight by approximating the 
chemical potential by its T = value. In this case, 
the variation of the critical temperature Tc with disor- 
der can be inferred from the gap equation through the 
relation p7ll28| 

where T^o = |Ao|e'*'/7r, with 7 ~ 0.577, is the critical tem- 
perature of the clean system, and ip{z) = T'{z)/T{z) is 
the di-gamma function. Again, considering explicitly the 



dependence of the bare order parameter | Aq | on the den- 
sity of excitations, it is possible to plot Tc/ TZy as a func- 
tion of both the density Pexflo ^^'^ disorder = 1/tc TZy 
(see Fig.O. The dependence of the critical temperature 
on the density at given values of the disorder closely mim- 
ics that of the order parameter |A|/7?.y given in Fig. 
Starting from values of the densities PcxO,Q ^ 1 the critical 
temperature starts to decrease as the Coulomb interac- 
tion weakens. When the density is high enough that the 
photon-like excitations are dominant, the critical temper- 
ature increases again and becomes insensitive to disorder. 

This concludes our discussion of the mean-field content 
of the theory. In the following section, in order to explore 
the spectrum of excitations above the condensate, we will 
develop an expansion of the low-energy action in fluctu- 
ations around the mean-field. 



IV. COLLECTIVE EXCITATIONS 

In order to address important physical characteristics 
of the condensed polariton system we need to go be- 
yond the mean field approximation and discuss excitation 
spectra. The excitations are visible through optical ab- 
sorption and photoluminescence, which can provide im- 
portant evidence for condensation in real systems. The 
phase stiffness of the condensate is another important 
parameter, which determines the transition temperature 
in the dilute limit - and in two-dimensional systems at 
all densities. In 2D, since there is no long range order 
(though there will be a Berezinskii-Kosterlitz-Thouless 
transition), the velocity of the Bogoliubov collective 
mode will determine the characteristic temperature scale 
at which a finite droplet of condensate can emit coherent 
light "2^. 

Thus in this section we will explore the impact of fiuc- 
tuations and how they translate to collective phenomena. 
The latter can be divided into separate contributions. In 
the first place, the system is susceptible to both spatial 
and temporal fiuctuations of the order parameters, jSj 
and \^\. In addition, there exist fiuctuations of the ma- 
trix fields Q(r) which encode the influence of mechanisms 
of quantum interference on the excitonic degrees of free- 
dom, the electrons and holes. In the following, we will 
consider an approximation engaging the rearrangement 
of the quasi-particle degrees of freedom in the condensate 
due to fluctuations of the order parameter at the saddle- 
point level. The impact of mesoscopic fiuctuations and 
weak localisation phenomena on the quasi-particle exci- 
tations is left as a subject for further investigation. Draw- 
ing on the analogy with the superconducting system, we 
note that related calculations of the pair susceptibility 
have been performed in Ref. [sOl in the context of the 
disordered superconducting film and, more directly, in 
Ref. [sJl in the context of a superconductor with mag- 
netic impurities. 

In section II Bl it was noted that the single particle 
Hamiltonian H1U|I can be expressed in terms of the sin- 
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gle complex composite order parameter A = S + gip. 
As such, in considering fluctuations of the total action, 
it is convenient to introduce a second independent field 
r = S — gip and express the fluctuations of the excitonic 
and photonic order parameters in terms of the fluctua- 
tions of A and F, viz. 



A(r,T) 
r(r,T) 



SA'-(r.T)+iSA^(r,T) , , 
' . !n (23) 



|A 

|r| + <5r^(r,T 



i<5r^(r, r) 



where |A| and |r| denote the mean-field expectation val- 
ues. Remembering that, at the mean-field level, the 
phases of S and ?/' are locked (see Eq. (^3): one can set 
|A| = iSKslV'l and |r| = while and are 

constrained by the relation ((TTjl . In Eq. longitudinal 
(l) and transverse (t) fluctuations represent respectively 
fluctuations of the amplitude and of the phase. Since our 
model is invariant under a global gauge transformation 
of the phase Hll|l . we can anticipate that the spectrum 
of fluctuations will exhibit one massless Goldstone mode 
and three massive ones. To develop a theory for fluctua- 
tions of the condensed field A, two of the three massive 
fiuctuations associated with the field F will be integrated 
out. 

Fluctuations of the order parameters A and F are 
not divorced from fluctuations of the fields Q. Indeed, 
through the coupling, one can expect the fluctuations of 
the order parameter to acquire a diffusive nature in the 
disordered system. Moreover, in the gapless phase, the 
existence of magnetic impurities should induce a mech- 
anism of dissipation. To capture these effects, one must 
determine the rearrangement of the field (or self-energy) 
Q due to fluctuations of the order parameter. 

Denoting the homogeneous saddle-point solution as- 
sociated with the mean- field as Q, let us represent the 
rearrangement due to fiuctuations as 



Q(r) = e^W/2Qe- 



W(r)/2 



(24) 



where the generators obey the symmetry relation 
{Q, Ty(r} = 0. Although the calculation of the rotation 
W in linear response is straightforward, the technical de- 
tails are involved and not illuminating. As such, they 
have been included as an appendix. Here, we note only 
that, at the level of the linear response, the longitudinal 
and transverse components of the order parameter fiuc- 
tuations do not get mixed by the Q-matrix action. By 
contrast, expanding S'ph + S's up to second order in (JA^-"^ 
and (5F'"''^, the (massive) fluctuations of the latter can 
be integrated out revealing a contribution which mixes 
the longitudinal and transverse response. Gathering all 
of the contributions together, the total action takes the 
form 

5tot = 5^'=[|A|,geff]+<5^tot , 

where S'*°[|A|,gcff] represents the mean- field Abrikosov 
and Gor'kov action for a superconductor with magnetic 



-4 



J 3 



G 0.5 




1.5 



\ y 




C = 0.2 






C = 0.5 






C = 0.8 






— C = i 






C = 1.5 








1 , 1 1 , 1 



H 3 



HO 

G 



0.5 





1 


1 




C = 0.2 






C = 0.5 






C = 0.8 






— C = i 


/ - 




C = 1.5 










1 , 1 




1 , 1 







FIG. 8: Zero temperature plot of the zero momentum longi- 
tudinal IIo^o ~ n|^),,o transverse nJ^Q — "^i^^^fi components 
versus the rescaled frequency Wf,/|A|, for different values of 
the disorder strength C = l/rc|A|. 



impurities HA5|I and bStot has the quadratic form 



Here the kernel G, , ^ „ denotes the inverse of the thermal 



Green's function or susceptibility: 



G 



H 



A, 



e 



A. 



-H^ 



(26) 



We refer to appendix ^ for details and the explicit ex- 
pressions of the Green's function elements. Let us note 
here that, if we restrict attention to the condensate phase 
|A| 0, the self-consistency equation H18|l translates to 
the condition 



eo,o 



(27) 



rendering the transverse fluctuations of the quasi-particle 
degrees of freedom massless. By contrast, when |A| = 0, 
the longitudinal and transverse kernels coincide, H'', „ = 
n|^^ p and acquire a mass, which, at zero temperature, is 
given by 



©0,0 - no;o = In 



IAoItc 
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As expected from the mean-field analysis, when 
2/|Ao|tc = 1 the mass disappears and the system un- 
dergoes a transition to a condensate phase. 

The form in Eq. I|25|) is positive definite, determin- 
ing the stability of the mean-field solution against fluc- 
tuations. Moreover it is already clear from the form of 
the kernel (|2(j|l and from the constraint (|27|l that the 
excitation spectra of the system admits a massless Gold- 
stone mode and a massive one. In fact, considering the 
zero frequency and zero momentum component, the anti- 
diagonal term Ag q vanishes together with the transverse 
one, {Gqq)227 while the longitudinal component (Gq q)ii 
remains finite. 



In fact, the precise functional forms IIJ;^'^p coincide 
with those found for a disordered superconductor in the 
presence of magnetic impurities 31]. Their frequency 
dependence is plotted in figure |S1 With H^j^^ q , no qual- 
itative change is associated with the transition to the 
gapless region for C > 1; the quadratic character of the 
dispersion at small frequencies is preserved. By contrast, 
as the system enters the gapless phase C ^ 1: with the 
longitudinal component n[^^ g exhibits a dissipative de- 
pendence scaling linearly with frequency. 

In order to elucidate this behaviour further, it is con- 
venient to affect a gradient expansion of the kernel H26|l , 
whereupon, at zero temperature, one finds 



g: 



' 1 I?p2 

^^ + 2"^W 



1 \^h\ 



2C |A| 



lAI 



1 DpP- 
7:a2-i-rr 



(28) 



with 



fli = %(() + a. 



a2 = ariO + o. 
b2^]:+b . 



(29) 



The coefficients ci = Cl(C), aL,T(C) and 6r t(C) derive from 



the expansion of the kernels IIJ;}^ p j32| . Although the 
explicit expression for their dependency on the disorder 
strength C can be rigorously derived (or equivalently can 
be numerically inferred from figure l|SIl), we note here 
that, as soon as |A| ^ 0, these parameters are of 0(1) 
and in particular, at zero disorder strength ^ = 0, they 
are given by c^{Q) = 1, a^{Q) = 7r/4, aT(0) = 7r/2, 6^(0) = 
1/6. More significantly, in the gapless region C > 1, the 
longitudinal kernel acquires a dissipative term and the 
coefficient 6l(C) changes of sign and becomes negative. 
In contrast, making use of the dimensionless parameters 
introduced in Eq. ©, the expansion of the kernels Qui^.p 



and A,, 



gives 



b 
I 



7^y 
7^y 



[g 



-gc{x)x\- 

S2 



[g-^ + ~gc{x)xf 



the zeros of the expression det(G^^^ p) 



Let us consider 

first the gapped phase, C < 1: where the determinant of 
the kernel can be evaluated exactly: 



det(Gj,^,p) 



El 



OJh 



El 



The energies -Ei,2, which can be shown for each value of 
the momentum to be real, are plotted in figure]^ From 
this we can infer that the spectrum of excitations, which 
is obtained by analytically continuing to real energies 
(i.e., substituting tj/j i— > ie), is given by a linear (massless) 
branch and a quadratic (massive) branch, where, expand- 
ing for small values of momentum, the gap of the mas- 
sive mode and the 'velocity' associated with the massless 
mode are respectively given by: 



Ei_ 2b2Ci+Ap 
lAI p^oV 



bib2 



a2Ci 




E2 



|A| V 62C1 + 2/2 



In order to gain some physical insight into these pa- 
rameters, once they are expressed in terms of the coef- 
ficients ()29|l. two opposite limits can be identified. In 
the first case, when 1 < poxOg ^ d/Air (i.e. for x » 1, 
gc{x)x ^ g and |A|/7?.y < 0(1)) one expects Ol.t ^ a 
and 6l.t S> fc, /, whereupon 



A. Gapped Phase 

The spectrum of the collective excitations is deter- 
mined by the poles of the matrix Gui^.p and therefore by 



E, ,/-±|A| 2V3|Ag| 



E2 ^ ./-^v^IpI - ^/^v/^|p| 
p^o y Ot C— »o 
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FIG. 9: Dispersion of the collective excitation spectrum in 
the gapped C < 1 (left) and gapless C > 1 (right) regions. 
Here, we have chosen the particular values ci = 0.8, ai = 7.1, 
a2 = 7.7, bi — 1, b2 — 1.2 and / — 0.6 (corresponding to 
pexOo ~ ^Q, £, ~ ^/'''cTZ-y ~ 1 and Dujc/t? ~ 0.1) in the first 
case, and ci = 0.1, ai = 7.6, a2 = 7.6, 6i = 0.4, 62 = 1.1 and 
/ = 0.7 (corresponding to pexOo = 16, ^ = 4 and Duc/c^ ~ 
0.1) in the second (note that this choice masks the actual 
linear dispersion of the massless mode). When ^ > 1, both 
the massive and massless excitations acquire a finite life time 
with, respectively the first finite and almost constant, while 
the second exhibits a quadratic dispersion ajid therefore null 
at zero momentum. Their line widths, Ki5i,2 ± 9-Ei,2, are 
explicitly plotted. 



Therefore, in the gapped phase, for values of the density 
at which the excitations are largely electronic in character 
(see Fig. 121), the gap for the massive collective excitations 
is proportional to the composite order parameter |A|, 
while the velocity of the massless excitations is give n by 
the product of the coherence length ^coh = \/£'/|A| and 
the order parameter [s^. 

In the opposite limit, PexOo ^ d/Ai: (i.e. for a; 0, 
and, taking into account the density dependence of the 
order parameter in this limit H22|) . for d/A-K > g'^/gdx)), 
we can use the approximations a ^ Ol.t and b, f ^ 6l,t, 
which translates to 



^1 ~ — r^/cy 

P-^O gc(X) 

When the density is so high that photons dominate, but 
still gc{x) 7^ 0, the gap in the collective excitation spec- 
trum is given by the ratio of the square of the photon 
coupling to the Coulomb coupling constant, while the 
velocity is proportional to the velocity of light c. Note 
that in absence of Coulomb interaction, gc = 0, the coef- 
ficient b vanishes and the situation is different: the gap 
is in fact proportional to |A| = gji/'l in the first limit and 
to IV'P/i'i^ in the very high density regime. 



B. Gapless Phase 

Let us consider now the gapless region C > 1- In the 
leading approximation, one may neglect the quadratic 
term 6iw^/|Ap, obtaining 

|A| pSo V^2Ci + 2/2\/ |A| +'4C(62Ci-l-2/2) |A| ' 

Therefore, analytically continuing to real energies 
(|w/i| ie), the massless modes E2 acquires a positive 
imaginary part with a quadratic dispersion in terms of 
the momentum p. In other words the massless mode has 
a finite life time (or line width) whenever p ^ 0. In or- 
der to evaluate the life time of the massive modes Ei, 
one has to restore the quadratic term feiw^/IAj. Consid- 
ering for simplicity a value of the density PexQo ^'^'^ 
disorder ^ = 1 /tc TZy for which the coefficient bi is over- 
all positive, the spectrum is plotted in figure |5| (Note 
that close to the transition between the gapped and the 
gapless region, the coefficient &i can become negative, in 
which case, to reestablish stability, one has to introduce 
the next order in frequency of the expansion.) 

We finally note that, in the non-condensed region 
|A| = (i.e., at zero temperature, for values of the dis- 
order strength 2/tc|Ao| > 1), the kernel acquires a mass 
proportional to ln(2/Tc|Ao|), while the massless modes 
disappear. 

C. Number of Photonic Occupied Modes 

Experimentally, the spectrum of collective excitations 
can be accessed indirectly by measuring the angle re- 
solved photoluminescencc intensity associated with the 
small leakage of photons from the cavity. The latter pro- 
vides a convenient way of mapping the occupation den- 
sity of the photon modes, 

AAph(p) = lim i^liT + t])iPp{t)) , 

where the in-plane quantum well momentum is deter- 
mined by the angle cf> at which the light is detected, 
|p| = (wc/c) tan^. Taken at the level of the mean-field 
approximation, the transition to the condensed phase 
would be associated with the development of a sharp co- 
herent peak in the density distribution at p = 0. How- 
ever, in the quasi two-dimensional quantum well geom- 
etry, the proliferation of massless phase modes destroys 
true long range order. Instead, the bulk transition to 
the condensed phase is of Berezinskii-Kosterlitz-Thouless 
type |3^. In this case, one expects the bulk transition to 
be reflected in the development of power-law correlations 
in the density distribution. Strictly speaking, to infer the 
characteristic power-law dependence of the correlations. 
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one should restore gauge invariance of the action under 
global phase rotations. However, guided by a parallel in- 
vestigation of trapped excitons in Ref . , one can infer 
the nature of the long ranged correlations simply from 
the spectrum of the phase mode determined above. 

In the present case the calculation is complicated 
slightly by the fact that, to establish the impact of the 
massless collective fluctuations on the photon distribu- 
tion, one must find their projection onto the phase mode 
associated with the photon field. This, in turn, intro- 
duces a 'form factor' which reflects the superfluid stiff- 
ness. Then, when parameterized as ip{r) = |V'(r)|e'^^'"-', 
discarding massive fluctuations, one obtains the charac- 
teristic power-law divergence |33 



p^op^(pL)2 

with the (superfluid) stiffness given by 



(30) 



quantum well gives rise to different domains in the phase 
diagram. Changing the density, we have shown that the 
character of the condensate changes from being excitonic 
in character to photonic. The effective coupling constant 
is determined by the sum of the direct screened Coulomb 
interaction and the effective interaction induced by the 
exchange of photons. When the condensate has largely 
excitonic character, we have shown that the 'charged' dis- 
order presents a pair-breaking channel which destabilises 
the condensate phase. Increasing the density of excita- 
tions, the system moves towards a region dominated by 
photons. Here, the increasing insensitivity of the con- 
densate to disorder inhibits the complete quenching of 
coherence, while leaving open the possibility of a gapless 
region where the condensate manifests the properties of a 
semiconductor laser. In the condensed phase, across the 
whole range of densities, we have determined the spec- 
trum of collective excitations. 



— — y > ^y'^ 



where fcoh = \/D/\A\ is the superconducting coherence 
length. The corresponding spatial correlation function 
assumes the characteristic logarithmic divergence at long 
distances: 



0{o)r 



where oc ^coh|A|/A:Br represents the thermal length. 
Restoring gauge invariance, and noting that the dom- 
inant contribution to the amplitude mode comes from 
mean-field, we obtain: 



A/'ph(p) 



ip-r-g(r) 



This integral can be exactly evaluated in two-dimensions, 
giving, for temperatures rj — k^T/2TTps < 2, the result: 



T(l-r7/2) 1 



r(r?/2) 



(31) 



As noted at the end of section HV Bl when |A| — 0, the 
amplitude and phase modes coincide and, in particular, 
the phase mode acquires a mass proportional to the dis- 
tance from the critical point. Here, in the non-condensed 
phase, the density distribution asymptotes to a constant 
at small p. 



V. CONCLUSIONS 

In this paper, we have investigated the properties of a 
model for cavity polaritons in the high density regime, 
where the interplay between Coulomb interaction, the 
photon mediated coupling and the effect of disorder in the 



The parameter regimes we have explored are in prin- 
ciple accessible in semiconductor microcavities. In III-V 
and II- VI semiconductors, microcavities have been con- 
structed with measured values of the Rabi splitting ap- 
proaching 30meV j^Land even larger values for an or- 
ganic semiconductor [3. Since the exciton Rydberg and 
inverse density of states take values similar to or smaller 
than this energy (depending on material), the dimension- 
less photon coupling g can be made of order unity. To 
reach the high density 'plasma' regime for the electron- 
hole system, one must obtain electron-hole densities PeiOo 
of order unity or greater, but this is indeed the typical op- 
erating density of a semiconductor laser. Indeed semicon- 
ductor lasers are conventionally modelled as a degenerate 
electron-hole plasma, though of course their operation is 
conventionally at temperatures somewhat larger than the 
Rydberg. Perhaps more important than their tempera- 
ture is that real lasers are open systems and thus have 
dynamical decoherence and dephasing processes that we 
have not included here. To the extent that pair-breaking 
by static disorder mimics dynamically induced decoher- 
ence, we may interpret the pair-breaking parameter ^ as 
the dephasing rate Tj"^ introduced into laser rate equa- 
tions. With this identification parameters for a conven- 
tional laser are indeed well in the limit C ^ 1- 



As a final remark, for the clean system, the model con- 
sidered in this work bears similarity with that discussed 
in the context of atomic condensation. The polariton 
Hamiltonian considered in this work has been used by 
many to model the atom/molecule interaction in a Fermi 
system positioned close to a Feshbach resonance (see, 
e.g., Ref. H^). With an appropriate dictionary of cor- 
respondence, both the mean-field content of the theory 
as well as the spectrum of collective excitations can be 
compared in the two systems. 
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APPENDIX A: SECOND VARIATION 

As outlined in section llVl the second variation 6Stot — 
SS[Q]+SSpYi+s of the total action is given by two contri- 
butions, one coming from the action for the semiclassical 
Green function (|13|l and the other from the photonic ^ 
and excitonic © actions. Starting with the first con- 
tribution, we consider the parameterization l|24(l for the 
Q{r) matrix and linearize the equation of motion (|14|l in 
the fluctuation terms W{r) and SA^^^^lr), 

- DV^W + [O, QW] - -^[Q, asa^'QWasaf'] 

ZTc 

= -[(5AV2 + 5AVi,Q] , 

where the operator 

6 = ia^af - |A|ct2 + ^^-iofQa-iuf , 



generates the mean-field equation of motion, [0,(5] = 0. 
The equation for Wiv) can be easily solved in term of the 
components W = WQ^+y^-cruf '' . Then, this solution can 
be substituted back into the action (|13|l . which expanded 
up to the second order, reads as 5[Q] = S\Q\ + 
where 

55[Q] = -^ y drtr([<5AV2 + <5AVi,Q]VF) . 
Substituting the solutions for and w, one obtains 

+5A^-..,-pn^.,p'5A^.,p] , (Al) 
where the kernels p are given by: 



fi^'^(ui,M2) 



Vi + wfyiT^ 



(A2) 



Dp2 



U1U2 T 1 



Turning to the second contribution, one has to expand 
iSph + 'S's around the mean field value of the photon field 
1^1 and of excitonic order parameter As suggested 
in section HVl it is convenient to introduce a new order 
parameter r(r, t) and parameterise the fluctuations as in 
Eq. (|23() . Once expanded up to the second order, the vari- 
ables 5T^^^^ can be eliminated via a Gaussian integral. In 

this way, fi-ph + S's = /?(wc - Ai)IV'P + /3|S|V5c + .JS'ph+E, 
where 

+^^-u,^~pQu.,.v5K,.p + '5A-.„,-pA-c.„p<5A^„p 
+<5AL^^^_pA^,^p<5A::„^,p} , (A3) 

and the terms O^^^^p and A^^,^ p are respectively given by: 

Q ^ 1 gM^(p)-M] + gc[^(p)-M]V.gc^fe 

UL^ {g^ + g,[u:{^) - + gl^l 



{g^ ^ cjA^iv) - ^Af + glu^l ' 

(A4) 



Note that the the term Xjvl? is here necessary since, 
in the thermodynamic limit, the Coulomb and photon 
coupling constant have to be rescaled with the volume 
according to the definitions ®. 

Finally, the mean-field value of the total action 

^«^[|A|,5eff]=/?(^c-M)lV'P+/3|S|V5c + ^[Q] , (A5) 

coincides with the Abrikosov and Gor'kov theory for a 
superconductor with magnetic impurities (see, e.g., t2d|). 
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